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ABSTRACT 

Motivated by the claimed discovery of a very massive galaxy (HUDF-JD2; M ~ 
5 X 10^^ Mq) at extreme redshift (z = 6.5) within the Hubble Ultra Deep Field 
(HUDF) (Mobasher et al. 2005), we have completed a systematic search for com- 
parably massive galaxies with z > 4 among the 2688 galaxies in our Ks < 
23.5 {AB) catalogue within the CDFS/GOODS-South field. This search was con- 
ducted using redshift estimates based on the recently-completed, uniquely-deep 11- 
band {B, i, z, J, Ks, 3.6/LtTO, 4.5/xm, 5.8/iTO, 8.0/im) imaging in this 125 square ar- 
cmin field, ~ 25 times larger than the NICMOS HUDF. To ensure completeness, our 
approach places no special emphasis on the standard y-drop, i-drop or z-drop criteria 
commonly used to pre-select candidate high-redshift galaxies. 

Initial spectral fitting, based on published catalogue SExtractor photometry, led 
us to conclude that at least 2669 of the galaxies in our sample lie at z < 4. This list 
includes several galaxies for which redshifts z > 4 have been previously proposed. We 
carried out a detailed investigation of the 19 remaining z > 4 candidates, performing 
aperture photometry on all images, and including marginal detections and formal non- 
detections in the fitting process. This led to the rejection of a further 13 galaxies to 
lower redshift. Moreover, subjecting HUDF-JD2 to the same analysis, we find that it 
lies at z ~ 2.2, rather than the extreme redshift favoured by Mobasher et al. (2005). 

The 6 remaining candidates appear to be credible examples of galaxies in the 
redshift range z = 4 — 6, with plausible stellar ages. However, refitting with allowance 
for extreme values of extinction we find that, even for these objects, statistically ac- 
ceptable solutions can be found at z < 3. In fact only 2 galaxies retain formally 
preferred high-redshift solutions. Moreover, the recently-released Spitzer MIPS imag- 
ing in GOODS-South has revealed that 5 of our 6 final z > 4 candidates are detected 
at 24/im. This was also the case for HUDF-JD2 (Mobasher et al. 2005), and provides 
further circumstantial evidence in favour of the moderate-redshift dusty solutions. We 
conclude that there is no convincing evidence for any galaxy with M > 3 x 10"'^^ Mq 
and z > 4 within the 125 square-arcmin GOODS-South field. We briefly discuss the 
implications of this null result, and revised expectations for the much larger (0.8 sq. 
degree), and deeper near-infrared UKIDSS Ultra Deep Survey now underway with 
WFCAM on the UKIRT. 
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1 INTRODUCTION 

Several hundred convincing galaxy candidates have now 
been uncovered at z ~ 5 — 6.5 (e.g. Bunker & Stanway 
2004; Taniguchi et al. 2005; Ouchi et al. 2005; Shioya et 
al. 2005; Yan et al. 2005; Bouwens et al. 2006), with a few 
(somewhat less convincing) galaxy candidates even reported 
at z > 7 (Bouwens et al. 2004). The discovery of such ob- 
jects has been used to set interesting new constraints on the 
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cosmic history of star- formation density {e.g. Stark & El- 
hs 2005, Bouwens & Illingworth 2006), but has not, as yet, 
presented a serious challenge to current theories of galaxy 
formation. This is because the masses of essentially all of 
these objects are relatively modest (M ~ lO^^Mo), and the 
observed large numbers of such objects are consistent with 
the predictions of at least some current galaxy-formation 
models (e.g. Nagamine et al. 2006). 

By contrast, the discovery of even a small number of 
very massive galaxies at these extreme redshifts can present 
a stern challenge for both semi-analytic and hydrodynamic 
galaxy-formation models. Indeed, given the steep decline in 
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the predicted number density of high-mass halos at z > 4 
(e.g. Somorvillo 2004), the discovery of a significant num- 
ber density of very massive objects at such redshifts has 
the potential to provide an interesting test of the now well- 
estabhshed paradigm of hierarchical structure growth within 
ACDM. 

For this reason, two recent studies have generated a lot 
of interest. First, Eyles et al. (2005), in their detailed study 
of 3 spectroscopically confirmed Lyman-break galaxies at 
z ~ -5.5 — 6. reported that these objects already contained 
a substantial, evolved mass of stars, apparently formed at 
redshifts as high as 2 ~ 7.5 — 13.5. Second, Mobasher et 
al. (2005) presented apparently convincing evidence for the 
existence of an extremely massive galaxy (HUDF-JD2; M ~ 
5 X 10" Mq) lying within the NICMOS Hubble Uhra Deep 
Field (HUDF) at the extreme redshift oi z = 6.5. Although 
this is just one object, its discovery within such a very small- 
area survey such as the HUDF is undoubtedly surprising and 
has already generated considerable interest {e.g. Panagia et 
al. 2005). 

Motivated by the discovery of HUDF-JD2, we decided 
to revisit our existing redshift determinations for Ks-hand 
selected galaxies in GOODS-South (Caputi et al. 2004; 2005; 
2006) , and to conduct a systematic search for the existence of 
any galaxies in this sample at very high redshift (2 > 4). The 
key point here is that HUDF-JD2 is sufHciently massive that 
it is relatively bright in the near-infrared, with Ks = 23.9. 
Thus it should be possible to detect comparable objects (less 
than a factor of 1.5 more massive at 2 ~ 6.5) within our 
complete Ks < 23.5 sample of 2898 objects (gala^cies and 
stars) in the GOODS-South field, which covers a solid angle 
~ 25 times larger than that subtended by the NICMOS 
HUDF. It is also clear that, in this extreme mass domain, 
the discovery of even one object within the entire GOODS- 
South field would be extremely important. 

Given the apparent evidence for some moderately 
evolved stellar populations in high-redshift galaxies, and be- 
cause we wished to search for any galaxies over a relatively 
wide redshift range 4 < 2 < 8, we decided not to base our 
search on strict colour criteria. Such criteria are frequently 
adopted in the selection of high-redshift galaxies, in order 
to ensure minimal contamination from low-redshift inter- 
lopers (e.g. Bunker & Stanway 2004). However, especially 
at 2 = 4 — 5, such clean selection inevitably comes at the 
expense of completeness, and introduces an inevitable bias 
in favour of the youngest, and hence lowest massdight ratio 
galaxies. Instead, we updated the dataset ultilised by Ca- 
puti et al. (2006) with the addition of the recently released 
ISAAC H-hand and complete Spitzer 4-band IRAC data, 
and then derived new redshift estimates for all 2688 galax- 
ies in the Ks < 23.5 sample by fitting a range of single and 
double component spectrophotometric models to the full 11- 
band optical-infrared photometry. 

The full results of this process are described in Cirasuolo 
et al. (2006) in which we have utilised this new dataset to 
explore the cosmological evolution of the galaxy mass func- 
tion out to z — 4. In the current paper we use this work 
simply as a starting point for the careful study of the rela- 
tively small subsamplc of galaxies which were found to have 
even marginally convincing solutions at 2 > 4. 

The paper is structured as follows. In Section 
2 we briefly describe how we have utilised the up- 



dated public dataset within the GOODS-South field 

to create a revised, complete sample of 2688 galax- 
ies with Ks < 23.5 and full, aperture-matched 11- 
band HST/ACS-hVLT/ISAAC+Spitzer/IRAC photometry. 
In Section 3 we describe the redshift-estimation technique, 
and demonstrate that it yields robust solutions at the cor- 
rect redshifts for known, spectroscopically-confirmed galax- 
ies at 2 ~ 5 — 6. Then, in Section 4, we describe the results 
of applying this technique to isolate a maximal sample of 
32 potential 2 > 4 galaxies within our GOODS-South sam- 
ple, (using 2.8-arcsec diameter aperture photometry) which 
was then refined down to a subsamplc of 19 serious candi- 
dates using smaller aperture SExtractor (Bertin & Arnout 
1996) measurements. Section 5 then presents the results of a 
detailed investigation of the mult i- frequency data for these 
remaining 19 candidates, with manual aperture photome- 
try leading to the rejection of a further 13 objects to lower 
redshift. In this section we also show that, by subjecting 
HUDF-JD2 to the same treatment, we find that it does not 
lie at z — 6.5 (as concluded by Mobasher et al. 2005), but 
is in fact a dusty, evolved galaxy at 2 ~ 2.15. The final in- 
terrogation of our remaining six candidate 2 > 4 massive 
galaxies in presented in Section 6. In particular, we explore 
the robustness of the derived redshifts when the assumed 
optical extinction is allowed to fioat up to values as large as 
Av = 6. Finally, in Section 7 we briefly discuss the impli- 
cations of our results in the context of theoretical models of 
structure formation. 

All optical and infrared magnitudes are given in the 
AB system (Oke & Gunn 1983). Masses and ages have been 
calculated assuming a cosmological model with JIm = 0.3, 
= 0.7, and Ho = 70kms"^Mpc"^ 



2 DATA AND SAMPLE DEFINITION 

2.1 New public data 

The survey area available for this study is determined by 

the coverage of the VLT ISAAC //-band imaging of the 
GOODS-South field released in September 2005. These data 
cover 125 arcmin'^, for which there already exists complete 
ISAAC imaging in J, and Ks, along with complete HST 
ACS imaging in B435 (3 orbits), Veoe (2.5 orbits), ^775 (2.5 
orbits) and zsso (5.0 orbits). 

Spitzer IRAC imaging, in all 4 IRAC bands has now also 
become available for the whole of this field, with the 3.6 /xm 
and 4.5 fxm imaging considered by Caputi et al. (2006) now 
supplemented by the 5.8 nm and 8.0 /xm imaging. 

2.2 Final sample properties 

The parent sample consists of the 2898 objects in this field 
which have Ks < 23.5 within a 2.8-arcsec diameter aperture. 
At this magnitude the limit the sample should be essentially 
100% complete. Star-galaxy separation (SExtractor stellar- 
icity parameter > 0.8 in the ACS 2-band images) led to the 
rejection of 210 objects as stars, leaving a final complete 
sample of 2688 galaxies. 

Following several recent deep multi-object spectroscopy 
campaigns (Vanzclla et al. 2005; Doherty et al. 2005; 
Vanzella et al. 2006; Roche et al. 2006), 850 galEixies in our 
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Figure 1. Spectral fits and versus redshift z for the four galaxies in the HUDF and GOODS-South field with measured spectroscopic 
redshifts z > 5, and published 7-band photometry as tabulated by Eyles et al. (2005) and Stark et al. (2006). Details of the model fits, 
and a comparison of estimated and spectroscopic redshifts are given in Table 1. 



sample now possess reliable spectroscopic redshifts. A fur- 
ther 188 galaxies, with R < 22.6, can be allocated solid red- 
shifts from the COMBO-17 survey (Wolf et al. 2004). Thus, 
for 1650 galaxies we required to derive an estimated redshift 
based on the available, deep, 11-waveband photometry. 

2.3 Initial multi- wavelength photometry 

The measurement of a robust 11-waveband spectral energy 
distribution for each source in the catalogue is non-trivial 
because the data span such a wide range in resolution (from 
~ 0.1 arcsec in the HST _B-band images, through to ~ 2 arc- 
sec in the Spitzer IRAC 8/im imaging). To maximise sensi- 
tivity at short wavelengths, it is tempting to consider the use 
of small (< 1 arcsec) apertures for source photometry in the 
ACS images. However, while most of the flux from genuine 
high-redshift {z > 4) galaxies is likely to lie within a 1 arcsec 
aperture, there are obvious dangers in combining small aper- 



ture HST photometry with the larger apertures required to 
obtain a robust estimate of total flux in the near-infrared 
ground-based imaging, and the Spitzer IRAC imaging. In 
particular, unless sufficient care is taken with the aperture 
corrections, the strength of any putative break between the 
J and zs5o filters can be exaggerated, leading to an erroneous 
conclusion in favour of a very high-redshift galaxy. 

To minimize the effect of any such bias we based our 
initial photometric catalogue on the use of 2.8-arcsec diam- 
eter apertures for all the optical and near-infrared data. For 
the 4 Spitzer IRAC bands we also used a 2.8-arcsec diameter 
aperture but then applied an aperture correction to estimate 
the anticipated total flux for a point source (ranging from 
0.55 mag. at 3.6^m to 1.0 mag. at 8^m). As we explain 
below, once we had used these data to isolate the subset of 
potential z > 4 galaxies, we considered alternative strategies 
to push the imaging closer to its photometric limit. 
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Table 1. A test of a our redshift estimation code at high redshift. Best fitting model parameters are given for our spectral fits to the 
broad-band photometry of the four spectroscopically confirmed z = 5 — 6 galaxies in the HUDF/GOODS-South field for which 7-band 

photometry has been published by Eylcs ct al. (2005) and Stark ct al. (2006). For each object we give the best fitting parameter values 
for both the best-fitting model with an exponentially declining star-formation rate, and for the best-fitting simple burst model. The fits 
arc plotted in Figure 1, along with versus estimated redshift, marginalised over age, normalization, varied star-formation history, and 
extinction. 



Source 


Zspec 


Zest 


Model type 




Age / Gyr 


Av 


Mass / 10^° Mo 


SBM03#1 


5.83 


5.70 ±0.07 


r 


= 0.3 Gyr 


5.26 


0.64 


0.0 


3.2 






5.40 ±0.15 




Burst 


11.63 


0.10 


0.0 


1.4 


SBM03#3 


5.78 


5.70 ± 0.07 


T 


= 1.0 Gyr 


9.01 


0.91 


0.0 


6.5 






5.55 ±0.15 




Burst 


19.26 


0.06 


0.0 


1.6 


GLARE#3001 


5.79 


5.75 ±0.12 


r 


= 0.3 Gyr 


3.80 


0.06 


0.6 


1.1 






5.62 ±0.15 




Burst 


4.52 


0.06 


0.0 


0.3 


32.8020 


5.55 


5.26 ±0.10 


T 


= 0.3 Gyr 


9.64 


0.91 


0.0 


15.8 






4.76 ±0.15 




Burst 


36.23 


0.06 


0.6 


5.6 



3 REDSHIFT ESTIMATION 



3.2 Tests 



3.1 Technique 

The photometric redshift for each galaxy was computed by 
fitting the 11 photometric data points (from the i3-band 
to Sum) with synthetic galaxy templates. These templates 
were produced using the stellar population synthesis models 
of Bruzual & Chariot (2003), assuming a Salpeter initial 
mass function (IMF) with a lower and upper mass cutoff at 
0.1 and 100 Mq respectively (we also explored the use of 
the Jimenez et al. (2004) models, and the Chabrier (2003) 
IMF, but these yielded inferior fits to the data). 

A range of templates was constructed based on different 
assumed star-formation histories. Specifically, we considered 
i) a single instantaneous starburst with passive evolution 
thereafter, ii) exponentially declining star-formation rates 
with e-folding times in the range 0.3 < r{Gyr) < 15, and iii) 
two-component burst models (to cope with the possibility 
of more stochastic star- formation histories). In all Ccises we 
adopted solar metallicity. 

To account for the effects of dust reddening we adopted 
the obscuration law of Calzetti et al. (2000). Initially we 
allowed F-band extinction to range up to Ay = 2, but ulti- 
mately we explored the expanded range < Av < 10 (see 
section 6). We also added a prescription for the Lyman se- 
ries absorption due to the HI clouds in the inter galactic 
medium, following Madau (1995). 

Some additional information was also utilised to exclude 
unreasonable redshift solutions. First the photometric red- 
shift was constrained to lie at z < 2 if the source had been 
detected in the (relatively shallow) f7-band imaging of the 
CDFS undertaken as part of the ESO Wide Field Survey 
(Arnouts et al. 2001). Second, a high-redshift solution was 
excluded if it resulted in a galaxy lying more than 1.3 mag- 
nitudes brightward of the K — z relation defined by the most 
luminous radio galaxies ( Willott et al. 2003) , in effect equiv- 
alent to a galaxy more massive than a present-day lOL* el- 
liptical. 



Our redshift estimation code has been tested for the ~ 1000 
galaxies in our GOODS-South /sTs-selected sample which 
possess a reliable spectroscopic or COMBO-17 redshift. As 
detailed in Cirasuolo et al. (2006), the la uncertainty in es- 
timated redshift inferred from this test is 6z/{l + z) ~ 0.07. 

However, of more specific interest for the present study 

is the ability of our method to accurately estimate the red- 
shifts of known objects at 2 > 4. To check this, we applied 
our code to estimate the redshifts of the four spectroscopi- 
cally confirmed « ~ 5 — 6 galaxies in the UDF/GOODS field. 
This is a useful test because the published multi-wavelength 
photometry for these 4 sources (Eylcs et al. 2005, Stark et 
al. 2006) provides comparable wavelength coverage to the 
photometry available for our complete GOODS-South sam- 
ple. 

The results of this test are summarized in Table 1 , and 
illustrated in Figure 1 which shows, for each galaxy, a plot 
of versus z (marginalised over age, normalization, and 
Av) and a comparison of the best fitting model SED with 
the broad-band photometry. 

Clearly, the code does an excellent job of recovering 
the redshifts of those galaxies. However, this is arguably a 
rather easy test, since these arc, by nature of their selection, 
rather clear-cut examples of young Lyman-break galaxies. 
As discussed further below, and as illustrated by the exam- 
ple of HUDF-JD2 (Mobasher et al. 2005), the situation is 
inevitably more confused if one is dealing with a potentially 
more evolved stellar population, with a less blue SED long- 
ward of the Lyman break. Nevertheless, the results shown in 
Figure 1 do at least confirm that our multi-wavelength SED 
fitting technique can efficiently and unambiguously identify 
high-redshift galaxies without recourse to any pre-selection 
of candidates based on, for example, specific colour criteria. 

For completeness, we also note that the masses and ages 
we have derived for these galaxies are in excellent agreement 
with those derived by Eyles et al. (2005) and Stark et al. 
(2006). 
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Table 2. Estimated redshifts and other model parameter values for fits to the observed spectral energy distribution of HUDF-JD2 using 
i) the published 11-waveband photometry (Mobasher et al. 2005) and ii) the revised photometry described in Section 5.2. 



Source Zest Model type Age / Gyr Ay Mass / lO" Mq 



HUDF-JD2: published photometry 6.20 ± 0.20 Burst 4.4 0.57 0.0 6.0 

HUDF-JD2; revised photometry 2.15 it 0.30 Burst 10.9 0.10 3.8 0.8 
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Figure 2. Progressive refinement of the estimated redshift distri- 
bution of the initial 32-source 2 > 4 candidate list as the photo- 
metric measurements axe pushed closer to the limit allowed by the 
imaging data. The top panel shows the redshift distribution based 
on initial 2.8-aj:csec diameter aperture magnitudes, and adopting 
formal limits for the non-detections in the ACS optical bands. 
The middle panel shows the impact of moving to smaller (1-arcsec 
diameter) SExtractor magnitudes. The bottom panel shows the 
effect of refitting the 19 objects which lie at z > 4 in the middle 
panel to the manual aperture photometry tablulated in Tables 
3 and 1. The 6 objects which apparently remain at 2 > 4 are 
discussed in detail in Section 5.3 and Section 6. 



4 REFINEMENT OF A HIGH-REDSHIFT 
SAMPLE 

4.1 Initial 2 > 4 sample 

Initial application of our redshift estimation code to the 
full 2688-galaxy 11-wavcband photometric datasot described 
above yielded formally acceptable solutions at 2 > 4 for 32 
galax;ies. This subset did not contain any of the objects for 
which redshifts have been obtained with the VLT (Vanzella 
et al. 2005, 2006), nor did it include any of the 18 claimed 
V-drop or J-drop galaxies within our sample listed by Bre- 
mer et al. (2004) and Dickinson et al. (2003), all of which 
were found to lie at 2 < 3. It also contains only 1 of the 7 
galaxies within our sample which are listed as having 2 > 4 
in the GOODS-MUSIC catalogue (Grazian et al. 2006). 

The initial derived redshift distribution for these 32 ob- 
jects is shown in the top panel of Figure 2. 

4.2 Revised 2 > 4 sample after l-arcsec catalogue 
search 

We next explored the 1-arcsec diameter aperture HST SEx- 
tractor catalogue for these 32 sources. This was done to 
check if any non-detections in the larger aperture became 
detections with reduced noise. The result of this was the 
relegation of a further 13 galaxies to lower redshifts, leaving 
a reduced subset of 19 potential 2 > 4 galaxies warranting 
further investigation. The result of this stage in the filtering 
is shown by the central panel in Figure 2, and details of the 
remaining 19 sources are given in Table 2. 



5 DETAILED STUDY OF THE FINAL 
HIGH-REDSHIFT SAMPLE 

5.1 Direct aperture photometry 

At this point in the analysis, the values adopted for limiting 
magnitudes in the case of SExtractor non detections become 

crucial. In particular, the extent to which a high-redsliift 
solution is favoured can be critically dependent on how any 
apparent non-detections are treated. 

One way to tackle this is the approach taken by 
Mobasher et al. (2005) in their analysis of HUDF-JD2. 
Faced with apparent non detection of this object in the deep 
HUDF ACS optical imaging, they adopted 2a upper limits 
of B = 30.6, V = 31.0, 1775 = 30.9, 2350 = 30.3, within their 
chosen 0.9-arcsec diameter aperture. We note here that at 
least some of these limits seem rather deep compared to the 
limits adopted by some other authors working with HUDF 
data. For example Bouwens et al. (2004), in their search for 
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zgso-drop galaxies in the HUDF, appear to have found a typ- 
ical 2a limit of zsso > 29.8 within a 0.6-arcsec diameter aper- 
ture, which converts to zgso > 29.4 for a 0.9-arcsec diameter 
aperture, almost a magnitude shallower than adopted by 
Mobasher et al. (2005). 

We investigated the effect of calculating equivalent lim- 
its for our GOODS non-detections (i.e. scaling the Mobasher 
et al. values to a 1-arcsec diameter aperture, and correct- 
ing for the GOODS:HUDF orbit ratio in each band). This 
leads to the conclusion that our adopted 2a limits should be 
B = 28.9, V = 29.2, 1775 = 28.6, zsso = 28.3. 

However, as illustrated in the central panel of Figure 2, 
the adoption of such limits would lead us to the conclusion 
that many of our remaining 19 sources lie at very high red- 
shifts. Faced with such a radical conclusion we decided to 
abandon the inferred detection limits, and perform manual 
aperture photometry for all 19 galaxies to establish the true 
level of signal and noise, so that even the apparent SExtrac- 
tor non-detections could be properly incorporated within the 
fitting in a consistent manner. 

For consistency, we therefore decided to also subject 
HUDF-JD2 to the same type of analysis. 



5.2 HUDF2 revisited 

In the top panel of Figure 3 we shown our own, independent 
fit to the photometry for HUDF-JD2 published by Mobasher 
et al. (2005). Our preferred redshift of z — 6.2 ± 0.2 is con- 
sistent with the value of z = 6.5 derived by Mobasher et 
al. (2005), and we agree that the high-redshift solution is 
significantly favoured over the alternative option of a dust- 
obscured galaxy at z ~ 2 — 3. The parameter values for our 
best fit solution to the published photometry are given in 
Table 2. 

However, performing our own manual photometry on 
the public HUDF images, through a 0.9-arcsec diameter 
aperture (as adopted by Mobasher et al.) we derive B > 29.9 
(la), V = 29.6±0.35, 1775 = 29.5±0.5 and zgso = 28.6±0.35. 
None of these highly marginal detections could, in isolation, 
be described as very convincing. However, their cumulative 
effect on the best fitting solution is dramatic, as illustrated 
in the lower panel of Figure 3. The best fitting result is now 
at z — 2.15, and the high-redshift solution is formally ex- 
cluded. The parameter values for our best fit solution to this 
revised photometry are also given in Table 2. We conclude, 
therefore, that HUDF-JD2 lies at z < 3, and not at z > 6. 



5.3 Final photometry and redshift distribution 

The positions, and final 11- waveband photometry for the 19 
remaining z > 4 candidates are given in Tables 3 and 4. Op- 
tical magnitudes are based on manual photometry through 
a 1-arcsec diameter aperture as described above, with errors 
given in magnitudes if the signahnoise ratio is greater than 
3, and as a fractional error in fiux if the 'detection' is less 
significant. To avoid bias, the J, H, Ks values have been re- 
derived also using a 1-arcsec diameter aperture, with a min- 
imal (effectively stellar) aperture correction of —0.5 magni- 
tudes applied to correct for missing ffux. 

The estimated redshifts resulting from fitting to the fi- 
nal photometry for these 19 galaxies are listed in Table 5, 
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Figure 3. Spectral fits and versus redshift 2 for HUDF-JD2. 
The top panel shows the result of applying our model fitting to 
the photometric data published by Mobasher et al. (2005). The 
lower panel shows the result of fitting to our own independent 
photometry, as discussed in Section 5.2. Best-fitting model pa- 
rameters for both sets of data are given in Table 2. 



with the resulting redshift distribution shown in the bottom 
panel of Figure 2. Clearly all 19 of these galaxies are interest- 
ing objects, with a sharp decline in flux-density shortward of 
the J-band, but now only 6 galaxies retain credible solutions 
at z > 4. Plots of versus redshift (marginalised over age, 
extinction and star-formation history) and the best-fitting 
model SEDs for these 6 galaxies are shown in Figure 4. 

For comparison, we have also listed in Table 5 the es- 
timated redshifts for these 19 galaxies recently released by 
Grazian et al. (2006) as part of the GOODS-MUSIC project. 
We have also applied our own spectral fitting technique to 
the GOODS-MUSIC photometry for these objects, to allow 
us to explore the extent to which any disagreements may 
depend on photometry or model fitting. Inspection of Ta- 
ble 5 shows that, for many of these galaxies, the agreement 
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Table 3. Positions and 1-arcsec diameter aperture HST ACS optical magnitudes for the 19 candidate z > 4 galaxies in our GOODS- 
South sample. Errors are given in magnitudes if the signahnoise ratio is greater than 3, and as a percentage error in flux density if the 
'detection' is less significant. 1-a limits are given when the detected flux density was zero or negative (see Section 5.3). 
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2476 


03 


32 


37, 


.86 


-27 


52 


01, 


.3 


> 28.8 


> 29.0 


26.8 ±0.25 


26.1 ±0.20 


2507 


03 


32 


19, 


.67 


-27 


46 


02. 


.0 


28.9 (120%) 


> 29.5 


> 28.6 


27.8 (65%) 


2600 


03 


32 


38, 


.34 


-27 


51 


01, 


.0 


28.0 (70%) 


27.8 ± 0.35 


27.2 ± 0.30 


27.0 (45%) 


2609 


03 


32 


56, 


.10 


-27 


52 


05, 


,0 


28.3 (60%) 


28.1 (50%) 


27.1 ± 0.18 


26.3 ±0.20 


2694 


03 


32 


21, 


.99 


-27 


51 


11, 


,9 


28.3 (70%) 


28.7 (70%) 


27.7 (60%) 


28.9 (200%) 


2869 


03 


32 


17, 


.99 


-27 


50 


52, 


,7 


28.1 (70%) 


27.25 ± 0.25 


26.1 ± 0.15 


26.1 ±0.30 


2895 


03 


32, 


.16, 


,83 


-27 


49 


07, 


,7 


29.7 (250%) 


28.6 (50%) 


27.6 (60%) 


27.2 (70%) 


2957 


03 


32 


37, 


.06 


-27 


44 


19, 


,1 


28.8 (180%) 


28.9 (70%) 


28.4 (130%) 


26.9 (55%) 


2958 


03 


32 


42, 


.08 


-27 


41 


41, 


.3 


> 29.1 


> 28.8 


> 28.3 


> 28.6 


3021 


03 


32 


27 


.14 


-27 


53 


11, 


.6 


> 28.7 


28.2 (75%) 


26.7 ±0.35 


26.8 (35%) 


3037 


03 


32 


28, 


.21 


-27 


51 


16. 


.2 


> 29.1 


27.3 d= 0.20 


26.6 ±0.15 


26.4 ±0.25 


3048 


03 


32 


19 


.57 


-27 


41 


39, 


.9 


> 28.2 


> 29.1 


27.0 (50%) 


26.8 ± 0.33 


3087 


03 


32 


54, 


.81 


-27 


51 


38, 


,8 


27.8 (40%) 


27.8 (40%) 


28.5 (50%) 


26.8 ±0.25 


3088 


03 


32 


39, 


.13 


-27 


51 


05, 


.0 


28.5 (120%) 


28.9 (65%) 


26.8 ±0.15 


27.9 (60%) 


3122 


03 


32 


21, 


.69 


-27 


42 


42, 


.3 


> 29.1 


27.9 (40%) 


28.7 (150%) 


27.9 (80%) 



Table 4. Infrared photometry for the 19 candidate z > 1 galaxies. ,/, H, Ks magnitudes have been rederived from the public VLT ISAAC 
imaging using 1-arcsec diameter apertures, and applying a point-source correction of —0.5 magnitudes. The Spitzer IRAC magnitudes 
have also been corrected to account for the point-spread function, and typical errors adopted for all objects at each waveband. 



Name 


J 




H 




Ks 




3.6/im 




4.5;um 




5.8;im 




S 


!.0/im 




1865 


24.33 ± 0. 


.13 


23.52 ±0, 


.06 


22.62 ±0. 


,05 


20.96 ± 0. 


,20 


20.72 ± 0. 


.20 


20.64 ± 0. 


.30 


21, 


.02 ± 0, 


.30 


2028 


23.86 ± 0. 


,07 


23.15 ±0, 


.07 


22.66 ±0. 


,04 


20.55 ±0. 


,20 


20.77 ± 0. 


,20 


20.47 ± 0. 


,30 


20 


.60 ± 


.30 


2336 


24.98 ± 0, 


17 


23.86 ± 0, 


,12 


22.96 ±0, 


,04 


21.29 ±0, 


,20 


20.81 ±0, 


,20 


20.34 ±0, 


,30 


20, 


.30 ±0, 


.30 


2351 


25.10 ± 0, 


,34 


23.64 ± 0, 


,13 


22.98 ±0, 


,07 


21.33 ±0, 


,20 


21.10 ±0, 


,20 


21.04 ±0, 


,30 


21, 


.26 ±0, 


.30 


2476 


23.99 ±0, 


,08 


23.30 ±0, 


,08 


22.96 ±0, 


,05 


22.14 ±0, 


,20 


21.92 ±0, 


,20 


21.96 ±0, 


,30 


22, 


.26 ±0, 


.30 


2507 


25.48 ±0, 


,26 


24.41 ±0, 


,16 


23.04 ±0, 


,04 


21.20 ±0, 


,20 


20.77 ±0, 


,20 


20.25 ±0, 


,30 


20, 


.43 ±0, 


.30 


2600 


25.20 ±0, 


,22 


23.83 ±0, 


,09 


22.93 ± 0, 


,05 


21.59 ±0, 


,20 


21.26 ±0, 


,20 


21.23 ±0, 


,30 


21, 


.48 ±0, 


.30 


2609 


24.63 ± 0. 


,14 


23.42 ± 0, 


.08 


23.12 ±0. 


,06 


21.78 ±0. 


,20 


21.47 ±0. 


,20 


21.64 ±0. 


,30 


21, 


.71 ±0, 


.30 


2694 


25.87 ±0. 


,41 


24.59 ±0. 


.22 


23.08 ± 0. 


.05 


21.31 ±0. 


.20 


20.88 ± 0. 


.20 


20.47 ± 0. 


.30 


20, 


.78 ± 0, 


.30 


2869 


24.89 ± 0. 


,20 


24.06 ± 0, 


.15 


23.18 ±0. 


,06 


21.57 ±0. 


,20 


21.37 ±0. 


,20 


21.33 ±0. 


,30 


20, 


.41 ± 0, 


.30 


2895 


25.01 ±0. 


,18 


24.13 ±0, 


.15 


22.99 ± 0. 


,04 


21.68 ±0. 


,20 


21.57 ±0. 


,20 


21.73 ± 0. 


,30 


21, 


.85 ± 0, 


.30 


2957 


25.80 ± 0. 


,33 


24.65 ± 0, 


.15 


23.05 ± 0. 


.04 


21.89 ±0. 


.20 


21.52 ±0. 


.20 


21.06 ±0. 


.30 


21, 


.32 ± 0, 


.30 


2958 


25.51 ±0, 


,33 


24.56 ±0, 


.23 


23.34 ± 0. 


,09 


21.34 ±0. 


,20 


20.95 ±0. 


,20 


20.48 ± 0. 


,30 


20 


.44 ± 0, 


.30 


3021 


25.27 ±0, 


19 


24.29 ±0, 


,15 


23.46 ± 0, 


,07 


21.77 ±0, 


,20 


21.49 ±0, 


,20 


21.09 ±0, 


,30 


21, 


.19 ±0, 


.30 


3037 


24.60 ± 0, 


12 


24.01 ±0, 


,12 


23.39 ±0, 


,07 


22.38 ±0, 


,20 


22.13 ±0, 


,20 


21.53 ±0, 


,30 


22, 


.18 ±0, 


.30 


3048 


25.60 ±0, 


,31 


24.26 ±0, 


,23 


23.49 ±0, 


,08 


21.69 ±0, 


,20 


21.29 ±0, 


,20 


20.77 ±0, 


,30 


20, 


.93 ±0, 


.30 


3087 


25.30 ±0. 


,29 


24.16 ±0, 


.17 


23.22 ± 0. 


,07 


21.83 ±0. 


,20 


21.55 ±0. 


,20 


21.01 ±0. 


,30 


21, 


.26 ± 0, 


.30 


3088 


24.60 ± 0. 


,12 


24.19 ±0, 


.15 


23.33 ± 0. 


.08 


21.92 ±0. 


,20 


21.64 ±0. 


.20 


21.62 ±0. 


.30 


22, 


.03 ± 0, 


.30 


3122 


27.21 ± 0. 


,16 


23.85 ±0, 


.14 


23.28 ± 0. 


,09 


21.81 ±0. 


,20 


21.54 ±0. 


,20 


21.20 ±0. 


,30 


21, 


.35 ± 0, 


.30 



between all three rcdshift estimates is excellent. However, 
there are also obvious differences. Specifically, only one of 
these objects has a redshift z > 4 in the GOODS-MUSIC 
list, and our own solution for this object with either pho- 
tometry set lies at 2 < 4. Before discussing further the likely 
explanation for this apparent disagreement at high redshift, 
we describe below the result of our final analysis of our own 
remaining subset of 6 potential 2 > 4 candidates. 



6 MASSIVE GALAXIES AT HIGH REDSHIFT? 

The model fits, and plots of versus redshift shown in Fig- 
ure 4 were derived with dust extinction limited to Av < 2. 
As detailed in Table 6, most of these 'high-redshift' solu- 
tions are statistically acceptable, and yield plausible val- 
ues for many of the model parameters {e.g. ages less than 
the age of the universe at the epoch of interest). However, 
while obviously not as surprising as the claimed discovery of 
HUDF-JD2, the existence of 6 galaxies at 2 > 4 with masses 
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Figure 4. Spectral fits and versus estimated redsliift z for tlie six galaxies in our GOODS-South sample which still have plausible 
solutions at z > A after the sample refinement process described in Section 5. Best fitting model parameter values are given in Table 6. 
For Ay < 2, the high-redshift solutions shown here are formally preferred. 
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Table 5. Re-estimated redshifts for the 19 candidate z > 4 galaxies in our GOODS-South sample. Column 2 gives the estimated redshift 
(and associated uncertainty) which we have derived by applying the technique described in Section 3.1, to the revised photometric data 

tablulated in Tables 3 and 4. The values of minimum for these fits are given in column 3. In columns 4 and 5 wc give the corresponding 
object IDs and estimated redshifts recently published by Grazian et al. (2006), as part of the GOODS-MUSIC project. In column 6 we 
give a third estimate of the redshift for each object, this time applying our own technique to the photometric data published by Grazian 
et al. (2006). Column 7 gives the values of minimum for these fits, many of which are very large (apparently due to problems in fitting 
some of the claimed {/-band detections in the GOODS-MUSIC catalogue). 



Source Zesti (lu range) MUSIC ID MUSIC Zest Zest2 

(lo- range) x^ 



1865 


5.02 


(4.87-5.17) 


8.82 


30093 


2.04 


2.00 


(1.95-2.10) 


72.03 


2028 


1.95 


(1.80-2.10) 


42.30 


30120 


1.97 


1.85 


(1.80-1.95) 


60.69 


2336 


6.22 


(6.07-6.45) 


12.73 


30199 


2.73 


3.65 


(3.25-3.85) 


2.77 


2351 


1.88 


(1.80-2.25) 


5.15 


30142 


1.84 


1.85 


(1.80-1.95) 


25.34 


2476 


1.65 


(1.50-1.80) 


6.02 


30160 


1.59 


1.70 


(1.50-1.80) 


9.53 


2507 


4.87 


(4.72-5.10) 


2.20 


30049 


2.21 


1.90 


(1.70-2.20) 


2.13 


2600 


2.18 


(1.95-2.33) 


1.39 


30106 


2.05 


2.15 


(2.05-2.25) 


78.35 


2609 


2.33 


(2.18-2.40) 


10.79 


30146 


2.06 


2.15 


(2.05-2.25) 


27.48 


2694 


2.93 


(2.55-3.08) 


5.88 


30114 


2.57 


3.45 


(3.00-3.90) 


6.60 


2869 


3.30 


(3.08-3.45) 


10.08 


30115 


1.88 


1.55 


(1.45-1.65) 


65.62 


2895 


1.73 


(1.58-1.88) 


4.42 


30123 


1.94 


2.00 


(1.75-2.45) 


4.24 


2957 


3.45 


(3.30-3.83) 


3.81 


30048 


4.66 


3.75 


(3.50-4.30) 


1.49 


2958 


6.07 


(4.95-6.30) 


1.49 


30009 


1.73 


1.95 


(1.35-2.35) 


0.71 


3021 


3.00 


(2.78-3.23) 


8.33 


30175 


2.44 


2.20 


(2.10-2.30) 


21.98 


3037 


2.25 


(2.10-2.40) 


9.39 


5177 


2.73 


1.75 


(1.65-1.89) 


45.31 


3048 


4.87 


(4.72-5.10) 


5.83 


30036 


2.62 


2.45 


(2.10-2.70) 


3.64 


3087 


2.18 


(1.88-2.33) 


12.44 


30145 


2.53 


1.95 


(1.75-2.10) 


17.55 


3088 


6.00 


(5.85-6.37) 


14.14 


30097 


1.91 


1.45 


(1.40-1.55) 


14.48 


3122 


2.63 


(2.25-2.85) 


8.98 


30032 


3.54 


2.65 


(2.45-2.85) 


25.74 



M > 3 X 10 Mq within a 125 square arcmin field is still 
unexpected. 

Therefore, as a final step in the analysis we refitted 
these 6 galaxies with extinction now allowed to range up 
to Ay ^ 10. The results of this process are illustrated in 
the contour plots shown in Figure 5. Acceptable, alternative 
low-redshift solutions are now found for all 6 galaxies and, as 
detailed in Table 6, are now formally preferred for all except 
two galaxies (2507 and 3048), and even in these two cases 
the lower-redshift solution is formally acceptable. 

It clearly remains possible that one or more of these 6 
galaxies lies at extreme redshift, but a number of factors mit- 
igate against this conclusion. First, given equally acceptable 
solutions at ~ 2 (with moderate mass but high Ay) and 
z ~ 5 (with high mass and low Ay) the balance of other 
probabilities clearly favours the low redshift option. Sec- 
ond, during the completion of this work, the Spitzer MIPS 
24fj,m catalogue for the GOODS-South field was released 
to the public (Dickinson et al., in preparation). This cata- 
logue contains 24/im detections for 5 out of the 6 galaxies 
listed in Table 6 (all except 3088). Clearly, such a high de- 
tection rate at mid-infrared wavelengths strongly supports 
the lower redshift dusty solutions for these objects (we note 
that HUDF-JD2 was also detected at 24/xm; Mobasher et 
al. 2005). Third, our own re-analysis of HUDF-JD2, and 
the discrepancy between the GOODS-MUSIC redshift es- 
timates and our own results for several of these putative 
high-redshift galaxies, serves to demonstrate just how sensi- 
tive any conclusion in favour of 2 > 4 can be to the treatment 
of marginal and non-detections in the optical wavebands. As 
demonstrated by Figure 1, this is a much less serious issue 



for young/blue high-redshift candidates for which an accept- 
able lower-redshift dust-obscured solution often does not ex- 
ist. However, the results presented here for redder objects 
serves to illustrate just how hard it will be to unambigu- 
ously identify older (and hence potentially more massive) 
objects at extreme redshifts on the basis of spectral fitting 
to even the most comprehensive photometric dataset. 



7 DISCUSSION 

To illustrate the implications of this failed search for high- 
mass galaxies at « > 4 we show, in Figure 6, a plot of the 
comoving number density of galaxies with mass M > 3 x 
10" Mo within the GOODS-South field, as a function of 
redshift. This limit simply corresponds to the lowest mass 
found for the high-redshift solutions to the 6 galaxies listed 
in Table 6. 

The data points for the redshift bins < z < 1, 
1<2:<2, 2<z<3, and 3 < 2; < 4 are derived from the 
analysis of Cirasuolo et al. (2006) (updated from Caputi et 
al. 2006), while the reference point at z = has been derived 
by the appropriate integration of the galaxy mass function 
provided by Cole et al. (2001) (assuming a Salpeter IMF). 
The result of the unsuccessful search at z > 4 described 
here is illustrated by the upper limit plotted at 2 = 5, which 
represents the comoving number density in the redshift bin 
4 < 2 < 6 if one of the 6 final high-redshift candidates 
listed in Table 6 really does lie within this high-redshift bin. 
Also plotted in this figure are the comoving number den- 
sity of virialized dark-matter halos above 3 different mass 
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Figure 5. Likelihood contours illustrating the location of acceptable fits on the Ay-z plane for the 6 objects in our GOODS-South 
sample with putative redshifts z > 4. Allowing for large values of extinction, it can be seen that, for all 6 objects, alternative solutions 
exist at z ~ 2. Contours are shown at intervals of Ax^ = 1, 4, and 9 above the minimum value of x^- The best-fitting values of the model 
parameters for both the high-redshift and low-redshift solutions are given in Table 6. 



thresholds, as derived from modified Press-Schechter theory 
(Percival, priv. comm.). 

Several features of this diagram arc worthy of com- 
ment. First, at no redshift does the number density of high- 
mass galaxies present a fundamental problem for ACDM; 
the number density of potential dark-matter halos (10-20 
times more massive than the stellar masses of the galaxies) 
clearly exceeds the inferred number density of the massive 
galaxies. 

Nevertheless, out to redshift z ~ 3, the number density 
of these most massive galaxies changes only slowly (com- 
pared to the dark matter curves) and hence, for this class of 
galaxies, the inferred dark-matter:stellar mass ratio appears 
to evolve from ~ 80 at x ~ to ~ 20 at z ~ 3.5. Such 
inferred values do not appear unreasonable (see, for exam- 
ple, Mandelbaum et al. 2006), and this apparent evolution 
of the dark-matter:stellar mass ratio can be viewed as yet 



another manifestation of 'downsizing' or apparently 'anti- 
hierarchical' galaxy formation (e.g. Heavens et al. 2004). 

However, the upper limit derived here indicates that 
this 'anti-heirarchical' behaviour does not persist beyond 
z ~ 4, and that, at higher redshift, the number density 
of massive galaxies drops off as rapidly (or possibly more 
rapidly) than the number density of potential host dark- 
matter halos. 

An improved measurement of the comoving number 
density of these rare high-mass galaxies at high-redshift 

clearly requires a substantially larger survey than the 
125 square arcmin covered by the GOODS-South multi- 
wavelength imaging. The first such survey of the necessary 
depth and area is now underway. This is the UKIDSS Ul- 
tra Deep Survey (UDS), which covers 0.8 square degrees, 
and is designed to ultimately reach a 7^-band 5-a detection 
limit of AT = 25 (AB). This survey will therefore cover 25 
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Table 6. Best fitting model model parameter values for the final 6 candidate high-redshift objects. The values for the physical parameters, 
and the value of minimum is given for both the putative high-redshift (low Ay) and the moderate-redshift (high Ay) solution in each 
case. The locations of these solutions on the Ay — z plane are illustrated in Figure 5. 



Source RA (J2000) Dec (J2000) Zest (Ict range) Model type Age/Gyr Ay Mass/lO" Mq 



1865 


03 32 12.87 


-27 46 40.9 


5.02 (4.87 


- 5.17) 


Burst 


8.82 


0.56 


0.0 


6.0 








1.75 (1.55 


- 1.95) 


Burst 


5.20 


0.13 


3.6 


0.9 


2336 


03 32 25.25 


-27 52 30.3 


6.22 (6.07 


- 6.45) 


Burst 


12.73 


0.56 


0.0 


16.6 








2.20 (1.85 


- 2.35) 


Burst 


4.80 


0.13 


3.6 


2.0 


2507 


03 32 19.67 


-27 46 02.0 


4.87 (4.72 


- 5.10) 


Burst 


2.20 


0.57 


0.8 


11.2 








1.85 (1.55 


- 2.30) 


r = 15 Gyr 


2.43 


0.20 


4.8 


0.8 


2958 


03 32 42.08 


-27 41 41.3 


6.07 (4.95 


- 6.30) 


Burst 


1.49 


0.63 


0.2 


15.8 








2.05 (1.75 


- 2.40) 


Burst 


3.03 


0.18 


4.6 


1.9 


3048 


03 32 19.57 


-27 41 39.9 


4.87 (4.72 


- 5.10) 


Burst 


5.83 


0.31 


1.2 


5.9 








2.45 (2.10 


- 2.80) 


T = 0.3 Gyr 


2.39 


0.07 


3.6 


0.4 


3088 


03 32 39.13 


-27 51 05.0 


6.00 (5.85 


- 6.37) 


Burst 


14.14 


0.40 


0.0 


2.8 








1.90 (1.65 


- 2.00) 


Burst 


12.46 


0.20 


2.8 


0.6 



times the area of GOODS-South, to a magnitude limit sub- 
stantially deeper than the ifs-limit of the GOODS-South 
survey analysed here. A failure to find any galaxies with 
M > 3 X 1O"M0 and « > 4 within the UKIDSS UDS would 
move the upper limit shown at « ~ 5 in Figure 6 down 
by over an order of magnitude. A first analysis, based on 
the UKIDSS early data release, has provided evidence for 
a number of moderately massive galaxies at z > 5, but as 
yet has not revealed any as massive as M > 3 x lO^^M© 
(McLure et al. 2006). 
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Figure 6. A comparison of our best estimate of the evolution 

of the comoving number density of galaxies with stellar mass 
M > 3 X 10"*^^ Mq with the predicted evolution of the comoving 
number density of virialized dark-matter halos above 3 different 
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the analysis of Cirasuolo et al. (2006) (updated from Caputi et 
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the appropriate integration of the galaxy mass function provided 
by Cole et al. (2001) (assuming a Salpeter IMF). The result of 
the unsuccessful search at 2 > 4 described here is illustrated by 
the upper limit plotted at z = 5, which represents the comoving 
number density in the redshift bin 4 < z < 6 if one of the 6 fi- 
nal high-redshift candidates listed in Table 6 really does lie within 
this high-redshift bin. The curves show the comoving number den- 
sity of virialized dark-matter halos with M > 5 x IO'^'^Mq (solid 
line), M > 1 X IO^^Mq (dashed line), and M > 2.5 x IO^Mq 
(dotted line). The last (most massive) of these curves was deliber- 
ately chosen to coincide with the data-point at 2 = 0, and implies 
a dark-matter:stellar mass ratio of ~ 80 for these most massive 
galaxies in the present-day universe. 
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